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Abstract 

The Singwi, Sjolander, Tosi, Land (SSTL) approach is generalized to study the 
spin-correlation effects in a one dimensional electron gas. It is shown that the SSTL 
approach is capable of generating results comparable to the more widely used STLS 
approach. 
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1 Introduction 



The Singwi, Tosi, Land and Sjolander (STLS) approach|jI|] is a powerful theoretical tool 
in going beyond the random phase approximation (RPA) in studying the short range cor- 
relation effects of an interacting electron gas. It was originally developed for the three 
dimensional (3D) electron gas. In the STLS approach the short range correlation effects 
are described by the local field correction G{q) for the density response function. 

The STLS approach is later applied to the two dimensional and the one dimen- 
sional electron gas (1D)P]-|0 with long range Coulomb or short range interactions. 

The Lobo, Singwi and Tosi (LST) approach]^ was originally developed to calculate 
spin correlation effects in the interacting 3D electron gas. Although the calculated spin 
susceptibility is not in agreement with experiment, LST approach is applied to the two 
dimensional [^]-[ 11 1 and one dimensional electron gas|12|,[l^] problems. 

In this paper, we study the validity of yet another attempt to go beyond RPA, the 
Singwi, Sjolander, Tosi and Land (SSTL) approach ||14|| which was originally developed for 
a 3D electron gas. It is proposed as an improvement over the STLS approach to better 
take into account the compressibihty rule. For a review of the STLS, SSTL and LST 
approaches please see reference [|^. It is curious that though it is not very widely used, the 
SSTL approach is not investigated before for a low dimensional electron gas problem. It is 
therefore of interest to study the SSTL approach, to investigate its range of validity and to 
compare its results with the very widely used STLS approach. This is done in our earlier 
work for density correlations in a ID electron gas||16||. 

In the present paper, we generalize the SSTL approach to study the spin correlation 
effects in a one dimensional interacting electron gas. The organization of the paper is 
as follows: In section II we present the SSTL approach. The results for the structural 
properties are given in section III. A discussion of our results and the preformance of the 
SSTL approach is presented in section IV. 



2 Formalism 

The ground state of the electron gas at high density is paramagnetic. In a paramagnetic 
system the magnetic moments of the constituents are randomly distributed and their mag- 
netic moment is averaged to zero. If we apply a weak external magnetic field to such a 
system, it will develop a paramagnetic spin moment. The response of the system to the 
field can be studied via its wave vector and frequency dependent paramagnetic suscepti- 
bility function. In noninteracting electron system, the Pauli paramagnetic susceptibility 
term is pronounced instead of the paramagnetic susceptibility term. On the other hand, in 
an interacting electron system we have short range Coulomb and exchange effects. Then, 
the Pauli value must be improved. In this study, we investigate these effects within the 
self-consistent calculation approximation (SSTL) as in the case of density-density correla- 
tions. Furthermore, it is assumed that our system is a Fermi liquid and the electrons are 
embedded in a uniform positive background so that the system is neutral. 

In the self-consistent calculation approximation the wave vector and frequency depen- 
dent density and spin-density response functions, x'^(g,c<j) and x'*(g, cu) respectively, are 
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and 
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where Xo(?i^) is the free electron polarizabihty, ip'"^ {ip"-) is the spin symmetric (spin anti- 
symmetric) potential, g is the Lande factor and /^b is the Bohr magneton. 

The system responds to the applied magnetic field through the free particle susceptibihty 
modified by a local field correction. The static structure factor S{q)a,nd the magnetic static 
structure factor S{q) are related to the dynamic response functions by the fiuctuation- 
dissipation theorem as follows: 

S{q) = — du;Im{x'{q,u;)}, (3) 
nn Jo ^ ' 

and 
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nng'^/jL^ Jo 

In the mean field approximation, the effective potentials ip^{q) and ip"'{q) are defined as 

r{q) = V{q)[l - G{q)], r{q) - V{q)I{q), (5) 

where V{q) is the one dimensional Fourier transform of the Coulomb potential, G{q) and 
I{q) are the static local field correction arising from the short range Coulomb correlation 
and exchange-correlation effects for the density and spin-density responses, respectively. 
In the SSTL approximation they are given in one dimension by 



and 



n J (27r) q V{q) e{q) 

where n is the one dimensional electron gas density. The Fermi wave vector kp is re- 
lated to the linear (ID) electron density as n = 2gykp/Ti. The g^ is the valley degeneracy 
and in this work we assume gi, = 1, which is the case for the quantum wire structures 
GaAl / AlxGai-x-^s. In one dimension, the dimensionless density parameter is defined as 
Ts — n/AkpaQ, where a% — So/e^m* is the effective Bohr radius with background dielectric 
function Eq and effective electron mass m*. The e{q) is the static dielectric function which 
is given in terms of the density response function x'^(g) as 

^^=l + V{q)x'{q). (8) 
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In the STLS approximation, the potential under the integral sign in Eq. (6) and in Eq. 

(7) is not screened by e{q). 

The random phase approximation (RPA) describes the dielectric properties of the elec- 
tron gas very successfully at high electron densities. In RPA the short range correlations 
are assumed to be absent, i.e, the local field correction G{q) — 0. As the density of the sys- 
tem is lowered, the exchange and the correlation effects become important. The Hubbard 
approach (HA) was developed to improve the RPA by introducing a local field correction 
taking into account the exchange hole around the electron. In HA, the local field correction 
for spin correlations, Iniq), is given by 



lH{q) = 



(9) 



2 V{q) 

The spin symmetric and spin antisymmetric pair correlation functions are related to 
the static structure factor and the magnetic structure factor by a Fourier transform in any 
dimension respectively. In one dimension 



1 f°° 

9(1^) = 1 + o / dq cos{qr) [S{q) - 1] , 

Z JO 



and 
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9{r) = / dq cos{qr) 
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(10) 
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These can be written in terms of parallel spin pair correlation function, gii{r), and antipar- 
allel spin pair correlation function gii{r) as 



2 bnW 
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(12) 



The spin dependent potentials may be obtained by combining the V'''(?) and ■0"(?) ^ 

V^Tt(?) = r{q) + r{q), ^tt(?) = r{q) - r{q)- (13) 

The static density and spin-density susceptibilities in the self consistent calculation 
approximation may be obtained easily from Eq. (1) and (2) as 



X'iq) = 



p{SF)kF 


q — 2kF 
q + 2kF 
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q — 2kF 
q + 2kF 



(14) 



and 



X'{q) 



g^nlp{eF)kF 



In 


q + 2kF 
q — 2kF 
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q + 2kF 
q — 2kF 



(15) 
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where piep) = 2m* /irkp is the density of states at the Fermi energy in one dimensional 
electron gas. Note that the Pauli spin susceptibility is XPauU = iA^p{^f)W^- 



3 Results and Discussion 

The results presented in this section are obtained by solving Eqs. (4), (7) and (8) self- 
consistently. In Fig. 1, the magnetic structure factor S{q) is presented for different 
values for a wire of effective radius h = 2a^. The appearance of a peak at g = 2A;i;' as the 



density decreases has also been reported in earlier studies [0]. It gets progressively difficult 
to achieve convergence in the self-consistent calculation for > 1.8. An instability sets 
in at this value of r^, as also observed in earlier STLS studies [jl2|,[|l^. This unfortunately 
means that we are unable to study different possible phases of the ID electron gas as one 
varies the density or equivalently the parameter. 

The magnetic structure factor S{q) obtained using different approaches is shown in Fig. 
2. It seems that the Hubbard approximation has the most pronounced peak at shown values 
of Tg and the wire radius b. The sharpness of the SSTL peak comes out to be the least. 

In contrast to S{q), the static structure factor S{q) is easier to obtain even for larger 
values of r^, as observed also in earlier work[|l^]. The S{q) for = 2, and different values 
of the effective wire radii are shown in Fig. 3. 

The static local field corrections arising from the short range Coulomb correlation and 
exchange-correlation effects for the density and spin-density responses are shown in Figs. 
4 and 5, respectively. 

The spin dependent pair correlation functions g^^{r) and are shown in Figs. 6 and 

7, respectively, for different values. It is seen that has very weak dependence. 

The spin symmetric and anti symmetric effective potentials ipi-\{q) and ipi[{q) in units 
of V{q) are shown in Figs. 8 and 9. They are also compared with the STLS results. The 
effective potentials have rather similar behaviour in both cases. 

The static density susceptibility function is compared with STLS result in Fig. 10. 
The static spin-density susceptibility, on the other hand, is compared with STLS and HA 
results in Fig. 11. It is seen that the results in all three approaches are qualitatively and 
quantitatively similar. The peak at g = 2kp is less pronounced in the SSTL approach which 
seems to be its dominant character. 

Collective excitations in the one dimensional electron gas can be studied as the poles 
of the density and spin-density response functions x'^'^il,^)- The dispersion of the para- 
magnon peak, ujp{q) is shown in Fig. 12 for = 1.5. For small q, ujp{q) shows a linear 



behaviour, as also observed before in STLS work||12 . 

In summary, we have studied the spin correlations in a one dimensional electron gas, 
and have shown that the SSTL approach is capable of giving results qualitatively similar 
to those obtained by using the STLS approach. This is, of course, not enough to establish 
the performance of the SSTL approach in ID. The question of compressibility inconsistency 
should be settled before one considers the SSTL approach as a viable alternative to other 
approaches. 
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Figure 1: The magnetic structure factor S'(g) for different Vg values at wire radius h — 2a^ 
in the SSTL approximation. 




Figure 2: The magnetic structure factor S{q) at Ts = 1 and h — 2a ^ in different approxi- 
mations. 
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Figure 3: The static structure factor S{q) at rg ~ 2 for different wire radii in the SSTL 
approximation. 




Figure 4: The density local field correction G{q) in different approximations for rg — 1 and 



b = 2al. 
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Figure 5: The spin-density local field correction I{q) in different approximations for = 1 
and b — 2a r. 




Figure 6: The spin dependent pair correlation function g'ii-(r) for b — 2ag and several 
values in the SSTL approximation. 
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Figure 7: The spin dependent pair correlation function gii{r) for b — 2a% and several 
values in the SSTL approximation. 




Figure 8: The spin dependent effective potential ip-[-i{q) in units of V{q) for rg — 1 and 
b = 2a*^ in the SSTL and STLS approximations. 
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Figure 9: The spin dependent effective potential i/^uil) units of V{q) for 
b — 2a% in the SSTL and STLS approximations. 
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Figure 10: The static density response x'^(g) in units of p{sf) for Vs — 1 and h — 2a^ in 
SSTL and STLS approximations. 
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Figure 11: The static spin response x'^iQ) i^i units of g^fj.%p{6F) for = 1 and b = 2a% in 
different approximations. 




Figure 12: The dispersion of the paramagnon peak LJp{q) in units of Ep for 
b = 2a*^ in different approximations. 
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